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Abstract 


Single-point,  jet-fluid  concentration  measurements  obtained  from  high  Schmidt 
number  (Sc  ~  1.9  x  10^)  turbulent  jets  permit  an  investigation  of  temporal  scalar 
power  spectra,  for  jet  Reynolds  numbers  in  the  range  of  1.25  <  Re  x  10““*  <  7.2. 
At  intermediate  scales,  we  find  a  spectrum  with  a  logarithmic  derivative  (slope) 
that  is  increasing  with  Reynolds  number,  in  absolute  value,  but  less  than  s/j  at  the 
highest  Reynolds  number  in  our  experiments.  At  the  smallest  scales,  our  spectra 
exhibit  no  power-law  behavior,  possessing  a  log-normal  region  over  a  range  of 
scales  exceeding  a  factor  of  40,  in  some  cases. 
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1.  Introduction 

We  report  here  on  an  experimental  investigation  of  temporal  scalar  power  spec¬ 
tra  of  round,  high  Schmidt  number,  turbulent  jets.  It  is  i:>art  of  a  larger  effort  to 
better  understand  mixing  in  turbulent  free  shear  flows,  including  jets  and  shear  lay¬ 
ers.  In  these  experiments,  we  examined  the  jet  fluid  concentration  (scalar)  power 
spectra  for  several  reasons.  Spectra  are  sensitive  diagnostics  of  the  flow,  providing 
information  over  a  wide  range  of  scales.  Historically,  they  have  been  the  object  of 
a  great  deal  of  attention,  partially  because  it  is  possible  to  extract  predictions  for 
spectral  slopes  from  various  turbulence  theories  and  models. 

Key  among  these  turbulence  theories  are  the  1041  paper  by  Kolmogorov.* 
with  implications  for  the  scalar  field  in  the  inertial  range  discussed  by  Corrsin.^ 
Oboukhov,^  and,  for  higher  wavenumbers,  the  theory  by  Batchelor.'*  See.  for  exam¬ 
ple,  discussions  in  Monin  and  Yaglom,^  as  well  as  in  the  recent  review  by  Gibson.® 
Both  the  Corrsin  and  the  Oboukhov  theories  yield  predictions  of  power-law  spectra 
and  of  the  spectral  power-law  logarithmic  derivative,  or  slope,  as  it  will  be  subse¬ 
quently  referred  to  in  this  paper.  Specifically,  the  Corrsin  and  Oboukhov  theories 
predict  a  scalar  spectrum  proportional  to  in  the  inertial  range,  as  did  the 

1941  Kolmogorov  theory  for  the  energy  spectrum. 

For  energy  spectra,  this  has  been  observed  experimentally  under  many  condi¬ 
tions  (c/.  compilation  of  data  by  Chapman’^).  The  situation  is  less  clear  concerning 
scalair  spectra,  with  departures  from  the  predicted  behavior  continuing  to  fuel  de¬ 
bate  about  details  and  refinements  of  the  theory. 

Batchelor"*  and  Batchelor  et  a/.®  recognized  that  for  Schmidt,  or  Prandtl.  num¬ 
bers  away  from  unity  there  exists  an  additional,  scalar-diffusion,  scale,  now  com¬ 
monly  referred  to  as  the  Batchelor  scale,  which  admits  a  change  in  the  scalar  spec¬ 
tral  behavior.  The  Batchelor  theory"*  predicted  that  the  scalar  power  spectrum  at 
high  Schmidt  numbers  would  display  a  dependence  beyond  the  Kolmogorov 
wavenumber,  i.e.,  a  .spectral  slope  of  —1.  Measurements  in  the  laboratory  [e.g.. 
Gib.son  and  Schwarz^)  and  the  ocean  (e.y..  Grant  et  al}^)  were  subsequently  re¬ 
ported  to  be  in  accord  with  this  prediction. 
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On  the  other  hand,  more  recent  measurements  by  Gargett‘*  in  the  ocean  were 
found  not  to  exhibit  a  h~^  spectral  range  (see.  however,  discussion  by  Gibson  in 
Refs.  12  and  6).  The  same  result  was  noted  in  passive  scalar  mixing  mea.surements  in 
shear  layers^'^  and  in  scalar  measurements  in  turbulent  jets.*'*’^'’  Specifically,  despite 
adequate  resolution  in  those  experiments,  no  range  was  found  at  high  spatial 
wavenumbers,  or.  to  be  exact,  temporal  frequencies,  adding  to  the  questions  about 
the  universality,  if  not  the  validity,^®  of  the  classical  predictions  at  high  Schmidt 
number. 

The  issue  of  spatial,  vs.  temporal,  spectra  should  be  recognized  here.  The 
classical  theories  .died  deal  with  spatial  spectra.  One  could  argue,  therefitre.  that 
comparisons  of  measurements  of  temporal  spectra  with  predictions  of  spatial  spectra 
cannot  be  made  directly.  Two  points,  however,  should  be  noted  in  response.  First, 
the  overwhelming  majority  of  experimentally  obtained  spectra  reiK>rted  in  accord 
with  the  theoretical  predictions  have,  in  fact,  been  temporal.  Second,  at  least  in 
the  case  of  developing  flows  that  are  not  (statistically)  spatially  homogeneous  over 
the  range  of  spectral  scales  of  interest,  the  notion  of  a  spatial  spectrum  and  the 
assumption  of  a  statistically  spatially  homogeneous  turbulent  field  is  questionable. 
Temporal  spectra,  derived  from  point  measurements,  do  not  have  to  contend  with 
this  issue. 


2.  Experiment 

The  experiments  investigated  the  scalar  (concentration)  field  of  round,  axisym- 
metric,  momentum-dominated,  turbulent  jets  issuing  from  contoured  nozzles  into  a 
large,  quiescent  discharge  tank.  The  measurements  were  performed  in  the  far  field, 
on  the  centerline  of  the  jet.  Details  of  the  experimental  apparatus  have  appeared 
previously,*'* so  only  a  brief  overview  will  be  presented  here. 

The  experimental  facility  consists  of  three  major  parts:  the  jet  plenum,  noz¬ 
zle,  and  delivery  system;  a  large  reservoir  that  acts  as  the  discharge  tank;  and 
the  diagnostics,  consisting  of  an  argon-ion  laser,  focusing  optics,  collection  optics, 
detector,  signal-processing  electronics,  and  the  subsequent  data-processing.  The 


4 


working  fluid  is  water,  and  the  scalar  is  a  laser  dye  (disodium  fluorescein)  which  is 
homogeneously  premixed  with  the  jet  plenum  fluid.  The  resulting  Schmidt  number. 


with  u  the  kinematic  viscosity  of  water  and  “D  ~  5.2  x  10  *’cm*/s  the  estimated 
aqueous  species  diffusivity  of  the  fluorescein  dye  (Ref.  IS.  p.  280),  is  1.9  x  lOb 


The  jet  flow  was  established  and  maintained  by  pressurizing  a  downward  ori¬ 
ented  jet  plenum  with  gas.  Both  sonically-metered  and  blow-down,  nearly  constant 
pressure,  gas  delivery  configurations  were  used.  The  internal  exit  diameter  of  the  jet 
nozzle  is  2.5mm  (0.1  in.).  The  rectangular  discharge  tank  is  squaj-e  in  cross-section, 
approximately  2  m  high  and  Im  on  edge.  The  tank  bottom  is  over  600  nozzle  di¬ 
ameters  downstream.  Large  glass  windows  on  all  four  sides  provided  optical  access 
(see  Ref.  15  for  details). 

The  illumination  .source  was  an  argon- ion  laser  (Coherent  Innova  90).  The 
particular  unit  was  custom  selected  for  its  low  AM  noise  figure  (~  — 95dB)  over 
the  frequency  range  of  interest  in  these  experiments.  It  was  operated  at  a  power 
of  3.5  W  in  the  light-regulation  mode.  The  beam  was  .spatially  filtered,  expanded, 
collimated,  and  subsequently  focused  to  a  small  waist  located  on  the  centerline  of 
the  jet.  A  low  dye  concentration  wais  used  in  the  jet  j)lenum  (  10“®M).  with 

correspondingly  substantially  lower  concentrations  at  the  measuring  stations.  A 
more  detailed  discussion  of  this  and  related  issues  may  be  found  in  Refs.  14  and  15- 
The  emitted  fluorescence  intensity  was  then  proportional  to  the  local  scalar  (dye) 
concentration  c(x,  t ),  that  was,  in  turn,  averaged  over  the  extent  of  the  measurement 
volume. 

The  fluorescence  emitted  from  the  measurement  volume  was  collected  through 
a  narrow  slit  spatial  filter.  The  besun  profile  and  the  slit  width  defined  a  small, 
spatially-averaging  volume,  roughly  spherical  in  shape  and  of  extent  (diameter) 
fa  ^  50 /im.  as  estimated  by  direct  observation  using  a  cathetometer.  We  will 
return  to  tliis  quantity  later. 


A  photomultiplier  tube  (RCA  86451  was  used  to  detect  the  fluorescence  emitted 
from  this  volume.  Its  output  signal  was  amplified  by  a  custom-designed  low-noise 
transimpedance  amplifier,  low-pass  filtered  using  a  third-order  Butterworth  filter, 
digitized,  sampled  with  some  margin  with  respect  to  the  Nycjuist  frec}uency,  and 
stored  for  subseciuent  processing. 

The  measurements  to  be  discussed  here  were  made  in  the  far-  field,  on  the  axis 
of  the  jet,  for  jet  Reynolds  numbers  in  the  range  of 

1.25  X  IQ-*  <  Re  =  —  <  7.2  x  10‘  .  (2) 

u 

where  ito  is  the  nozzle  exit  velocity,  d  the  nozzle  exit  diameter,  and  u  is  the  kinematic 
viscosity.  Data  were  also  recorded  at  both  lower  and  higher  Reynolds  numbers.  The 
lower  Reynolds  number  jets,  however,  behaved  substantially  differently,  by  any  of  a 
number  of  criteria,  and  were  not  accepted  as  representing  bona  fide  turbulence.*  *  In 
the  other  limit,  the  jet  at  the  higher  Reynolds  number  { Re  =  10.2  x  10* )  produced  a 
distinct  hissing  sound.  This  was  probably  generated  by  the  transient  dilatation  and 
subsequent  oscillations  of  small  air  bubbles  caused  by  the  rapid  reduction  in  pressure 
in  such  bubbles  as  they  exited  the  nozzle,  or  by  cavitation  in  the  jet  near-field  region, 
or  both  (note  that  the  plenum  gauge  pressure  is  quadratic  with  Reynolds  number). 
See  discussion  and  references  in  Ref.  19.  pp.  452-453.  and  Ref.  20.  pp.  205-207,  for 
example.  As  a  result,  the  Re  =  10.2  x  10*  jet  was  exposed  to  different  near-field 
conditions  and  will  not  be  included  in  the  discussion  below.  See  Ref.  15  for  a  further 
documentation  of  the  data. 

Finally,  constraints  dictated  by  resolution  and  statistical  convergence,  vis-d-vts 
total  number  of  large  scale  structures  captured  and  length-of-run  considerations, 
led  to  measurement  stations  in  the  range, 

100  <  ^  <  305  ,  (3) 

d 


where  x  is  the  distance  from  the  nozzle  exit. 


c 


3.  Scalar  power  spectrum  estimation 

The  fluorescence  signal  0{t),  representing  the  photon  flux  incident  on  the  pho¬ 
todetector,  is  a  linear  function  of  the  spatial  average  of  the  converted  local  jet-fluid 
concentration  c(x,t)  o\'er  the  measurement  volume.  It  produces  a  signal  that  can 
be  approximated  by  a  convolution  over  c{t)  =  c(xo,/),  the  jet  fluid  concentration 
at,  say,  the  center  of  the  measurement  volume,  Xq,  be., 

<P(0  -  /  Kit  -  t')c{t')dt'  =  }u{t)®c{t)  .  (4) 

J  — OO 

In  this  expression.  /?a(0  models  the  impulse  response  of  the  spatio-temporal  aver¬ 
aging  process,  z.e.,  the  temporal  signal  that  would  be  measured  if  a  spatial  delta 
function  of  dye  was  converted  through  the  measurement  volume  at  the  local  flow 
velocity.  See  Fig.  1. 


C(X,  0 


Fig.  1  Sketch  of  jet  fluid  concentration  field  c(x,  t),  convected  through  the  mea¬ 
surement  volume  of  extent 

The  fluorescence  output  along  with  fluctuations  contributed  by  the  small 
laser  intensity  fluctuations,  convected  residual  non-uniformities  in  the  jet  plenum 
dye  concentration,  photon  shot  noise,  electronic  noise  generated  by  the  signal- 
processing  chain,  etc.,  was  processed  by  the  Butterwoith  low-pass  filter  to  produce 
the  total  signal 


.s(0  =  hi^pit)  ®  d){t)  +  nit)  =  h(t)<S)cit)  +  n(t)  , 


(5) 
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that  was  digitized  and  stored.  In  this  expression.  h{t)  =  /iLP('^)  0  is  the  total 
system  transfer  function  and  t7(/)  is  the  total,  low-pass-filtered  noise. 

Assuming  that  the  total  system  noise  n{f)  can  be  modeled  as  uncorrelated 
with  the  local  dye  concentration  time  history  c{t),  the  spectrum  Ss(u-’)  of  s{t)  can 
be  expressed  in  terms  of  the  spectrum  of  Q{t)  and  the  spectrum  5n(«-’)  of  the 

(low-pass  filtered)  noise  n{t),  be., 


~  +  Sn(^)  ,  (6) 

where,  from  Eq.  4, 

5e,{u,’)  ~  Sc{^')  ,  (7) 

with  Ha,{u))  =  !FT {  hs,{t]  },  the  Fourier  transform  of  /)a(0-  This  allorvs  us  to  relate 
the  total  signal  spectrum,  Sii(u;),  to  the  desired  scalar  fluctuation  spectnrm.  Sdu;). 
of  c(f),  i.e. 

S^u)  ~  |^f(a;)|2  5,(u;)  -f  5„(u;)  ~  Se(-')  +  5„(^)  ,  (8) 

where  =  :FT{  hit)  }  =  Hip(u;) 

For  these  experiments,  the  knee  of  the  Buttervvorth  low-pass  filter  was  set  sub¬ 
stantially  higher  than  the  range  of  frequencies  contained  in  5p(u;).  Its  main  purpose 
was  to  bandlimit  the  noise  and  de-alias  the  digitized  measurements,  aJlowing  the 
noise-floor  to  be  determined,  as  will  be  illustrated  in  the  spectra  presented  below. 
This  is  the  reason  the  modulus  squared  of  i^Lpfo;),  the  transfer  function  of  the 
low-pass  filter,  can  be  ignored  in  Eqs.  7  and  8,  and  wherever  it  multiplies 
and  5c(u;). 

Figure  2  illustrates  these  relations  by  comparing  the  spectrum  S.,{uj)  of  the 
total  signal  s(t),  i.e.,  fluorescence  4>{t)  plus  noise  n{t),  with  So(u;),  the  spectnim  of 
the  fluorescence  signal  alone.  The  latter  was  calculated  by  subtracting  the  estimated 
noise  spectrum,  S„(u>),  from  the  total  spectrum  Recalling  Eq.  6.  we  have 


5a^(u.-)  ~  5,(u;)  -  Sr,{uj)  . 


(9) 
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The  noise  spectrnm  was  assumed  to  be  white,  as  was  ftJund  to  be  the  case  in  sepa¬ 
rate  measurements  of  this  quantity  (see  also  Ref.  21  for  examples).  Nevertheless,  the 
result  is  not  sensitive  to  the  assumed  shape  of  the  noise  sj>cctrum  at  low  frequencies 
where  dominates.  The  data  processed  to  produce  the  si)ectra  in  Fig.  2  were 

recorded  at  x/d  =  100,  for  Re  =  1.25  x  10^.  Note  the  high  dynamic  range  of  the 
total  signal  spectrum,  i. e..  the  (logarithmic)  difference  of  tla'  low-freciuency  power 
to  noise-floor  power.  Note  also  that  the  span  to  one-half  the  (scaled)  sampling  fre¬ 
quency  is  well  beyond  the  noise-cross-over  freejuency.  As  can  be  seeui,  the  frequency 
extent  of  the  noise  floor  was  substantial. 


logjotf  Tg) 


Fig.  2  Sample  spectrum  of  the  total  signal  (solid  line;  fluorescence  -I-  noise)  and 
estimated  fluorescence  spectrum  (dashed  line:  fluorescence),  derived  from 
measurements  at  x/d  —  100.  Re  =  1.25  x  10^.  Frequency  scaled  by  Tf,{x}. 
the  local  large- sc  ale- passage  time. 


The  spectra  in  Fig.  2,  and  throughout  this  paper,  are  normalized  by  c*.  the 
square  of  the  local  mean  value  of  c{t).  multiplied  by  the  local  large-scale-passag(' 
time,  rf,{x),  and  plotted  in  terms  of  the  circular  frequency.  /,  scaled  by  Tf,{x).  In 

these  coordinates,  their  integral  produces  the  normalized  variance,  t-.g., 

j'i 


c  7o 


r 

'  ScifrsY 

L 

—9 

r** 

dfTf, 


(10) 
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The  local  large-scale-passage  time,  is  given  by, 


r«{;r)  = 


Wci(-t) 


where 

6{x)  ~  0.41  a-  , 


(11a) 


(11b) 


is  the  local  outer  scale  of  the  flow,  here  identified  with  the  (measured)  mean  trans¬ 
verse  extent  (visual  width)  of  the  conical  region  enveloinng  the  jet-fluid  (Ref.  15. 
Appendix  D).  and  Uc\{x)  is  the  mean  centerline  velocity.  The  latter  was  estimated 
from  the  relation 


Uci(.r) 


(11c) 


Uj  X  —  Xj 

where  uj  is  the  jet  velocity  and  Xj  the  jet  (virtual)  origin,  recommended  by  Chen 
and  Rodi.^^  This  spectrum  and  frequency  scaling  was  found  to  produce  similarity 
with  respect  to  the  downstream  coordinate,  x/d,  in  the  analysis  of  scalar  spectra 
measured  in  gas-phase  jets.^^ 


The  spectrum  5c(w)  of  the  scalar  fluctuations  c(t)  can,  at  least  formally,  be 
estimated  by  solving  Eq.  7,  he., 


‘5'c  (<*'’) 


S4^) 


or  Eq.  8,  yielding  a  result  in  terms  of  experimentally  estimated  quantities,  he., 

5,(a;)  -  5n(w) 


(12a) 


Sc(uj) 


(12b) 


We  should  note,  at  this  point,  that  Eq.  4,  assuming  a  fixed  ha.it),  is  not  a  proper 
equation  for  two  reasons.  First,  the  fluid  velocity  converting  the  c(x,  t)-field  through 
the  measurement  volume  is  not  a  constant.  Second,  different  scalar  field  lagrangian 
trajectories  through  the  measurement  volume  sample  chords  of  different  sizes  (and 
transit  times)  through  it.  These  two  effects  could,  in  principle,  be  expressed  as  two 
additional  convolutions  over  the  local  velocity  field  and  scalar  paths  through  the 
measurement  volume.  See  Fig.  1. 
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For  frequencies,  however,  corresponding  to  spatial  scales  of  the  order  of,  or  less 
than,  the  extent  of  the  measurement  volume,  where  the  /fal'*-’)  transfer  function 
will  have  an  effect  on  the  measurement  of  the  scalar  spectrum,  fractional  fluctua¬ 
tions  in  the  convecting  velocity  are  considerably  smaller  than  fractional  fluctuations 
in  c(x,t),  at  least  for  this  high-Schmidt-number  fluid.  They  are  also,  largely,  uncor¬ 
related  with  them.  Fluid  paths  through  the  measurement  volume  are  also  uncorre¬ 
lated  with  the  passively- convected  scalar  field.  As  a  consequence,  in  estimating  the 
spectrum  i.u  the  frequency  neighborhood  of 
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and  above,  these  two  effects  do  not  contribute. 


(13) 


Alternatively,  the  m  )st  general  expression  of  the  linear  dependence  of  the  flu¬ 
orescence  signal,  <p{t),  on  the  local  jet-fluid  concentration.  c(t).  is  given  by  Eq.  7. 
While  Eq.  7  follows  from  Eq.  4,  the  converse  is  not  true.  Eejuation  4  is  more  re¬ 
strictive,  also  requiring  a  definite  phase  relation  between  c{t)  and  o(t).  Fortunately, 
these  (unknown)  phase  relations  do  not  enter  in  the  relation  between  the  corre¬ 
sponding  spectra.  We  may  conclude  that  Eqs.  7  and  12  represent  valid  relations  for 
the  corresponding  spectra,  even  as  Eq.  4  cannot  be  accepted  as  a  correct  description 
of  the  time- history  of  the  fluorescence  signal  0[t). 


The  estimation  of  the  scalar  spectrum  5c(a-')  in  the  frequency  range  influenced 
by  the  spatio-temporal  averaging  of  the  measurement  process  also  requires  knowl¬ 
edge  of  the  Ha.iuj)  transfer  function.  This  can  be  estimated,  in  turn,  by  noting  that 
it  is  dominated  bj’  a  pole  corresponding  to  the  transit  time  of  the  flow  through  the 
measurement  volume,  he.. 


with 


1 

1  4-  lU! 


(14a) 


(14b) 


Performing  two  different  experiments,  under  as  identical  flow  conditions  as  was 
feasible,  at  two  different  spatial  resolutions,  we  were  able  to  compare  the  spectra, 
for  p  =  1,2  , 


Ssp{^)  ~  I *5'cf‘*^)  +  5’„p(u.’)  , 


( 15a) 
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corresponding  to  two  different  locations  of  the  doiniinuit  pole.  at.  say.  ta  =  Tj  and 
fa  -  f2  >  t-e..  for  p  =  1,2, 

//p(-c’)  ~  - - - -  .  (15h; 

1  +  ?  a/  Tp 


If  the  scalar'  spectrum,  Sf.(u.’),  could  be  tissuined  to  he  identical  in  the  two  experi¬ 
ments  (the  corresponding  noise  floors  were  determined  sejnuaNdy  in  each  case),  we 
see  that  the  ratio  of  the  two  fluorescence  spectra. 


G{u:\  Ti  ,  "2  ) 


■would  be  given  by 


G{a,';ri,r2)  ~ 


1  +  {j,'  To) 
1  -f-  {kJ  T\  ) 


1  , 

(f2/fl  )■ 


for  u,.'  <C  l/f2  : 
for  ^  1/t\  . 


(IGa) 


(ICb) 


independently  of  the,  as  yet  unknown,  scalar  spectrum,  Sc(u.’). 


lOSloltTj) 


Fig.  3  Dotted  lines:  Fluorescence  spectra  estimated  from  measurements  a.t  x/d  = 
100  and  Re  =  1.25  x  10'*  at  two  spatitil  resolutions.  Solid  line:  Estimated 
concentration  spectrum  (Eq.  12)  at  i/d  =  100  and  Re  =  1.25  x  10“*.  Dashed 
line:  Estimated  concentration  spectrum  at  x/d  =  305  and  Re  =  1.2  x  10^. 
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Two  fluorescence  spectra.  —  5,,,{u,’)  —  S„j{j.')  and  —  S„j(w'). 

from  one  such  pair  of  experiinenrs.  at  x/J  =  100  and  Re  =  1.25  x  iOh  are  {)lotted 
in  Fig.  3  (dotted  line.s).  The  tluor<'Scence  spectrum  with  the  larger  high  frequency 
content  is  tiu'  one  plotted  in  Fig.  2.  Power  spectra  were  com]>uted  numerically 
using  a  power  spectral  density  estinuition  inethodolog}'  that  has  evolved  ov(U’  the 
past  ten  years,  or  so.  A  d<jcumcntati<ju  of  some  of  its  earlier  h'atures  can  l>e  fouiui 
in  Ref.  24.  In  processing  the  data  in  the  experiments  reported  here,  the  j)ower 
spectral  density  estimation  program  computes  spectra  of  data  files  by  nn^ans  of 
FFT  methods,  and  incorporates  Hanning  windowing,  contiguous  record  ov«uIapping. 
and  parabolic  detrending,  among  other  features.  Records  up  to  2*'  points  can  be 
accommodated.  For  spectra  known  to  be  smooth,  the  program  can  i)rovide  third- 
octave  (~  1/10  decade)  gaussian  filtering,  sampled  at  20  points  per  decade,  to 
produce  the  final  spectra.  Thi.s  feature  was  used  for  all  the  s]>ectra  plotted  in 
this  paper,  not  so  much  for  smoothing,  but  to  reduce  the  number  of  points  to  a 
manageable  level  for  plotting  purposes  (note  that  2“'  =  G5.53G). 


logjotf  T^) 

Fig.  4  Computed  ratio,  G(u;;r].T2),  of  fluorescence  spectra  at  r/d  —  100  (Eq.  16). 
Circles;  Re  =  1.25x10’.  Squares:  Re  =  2.55  x  10'*.  Solid  line:  Least  squares 
fit  for  T]  and  r?. 
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The  ratios  G(^-:ri,r2)  for  a  ptiir  of  spectra,  measured  at  xjd  =  100.  at  He 
1.25  X  lO’^  (circles)  and  a  pair  at  He  =  2.55  x  10*  (squares),  respectively,  are  j)l<jtted 
in  Fig.  4.  .As  can  b(“  seen,  the  ratios  of  spectra  measured  at  different  Reynolds 
numbers  are  very  nearly  tlie  Scune.  in  accord  with  the  analysis  outlined  a}j(n-e.  <*v<*ii 
though,  as  we  will  see  below,  tlu'  spectra  themselves  are  He-de])endent , 

The  ctirve  (solid  line)  in  Fig.  4  is  a  least-squares  fit  for  tj  and  to  the  Iowct 
Reynolds  number  data,  that  were  characterized  by  the  liighcr  signal-to-noise  ratio, 
in  the  frequency  range  1.0  <  /r/,  <  2.S .  The  lower  limit  of  the  fit  range  is  chosen 
so  as  to  exclude  (the  small)  run-to-run  variations  at  frequencies  well  l>eIow  those 
affected  by  the  spatial  averaging.  The  upper  frequency  limit  is  dictated  by  the  less 
than  unity  signal-to-noi.se  ratic*  at  higher  frequencies  yet  (c/.  Fig.  2).  The  values  for 
rj  and  r2  e.^tiinated  by  this  pivuedure  were 

Str  Ti/r^  ~  4.2  X  10“^  and  2;r  r2/r/.  ~  C.G  x  10“'^  .  (17) 

respectively.  This  corresponds  to  an  effective  spatial  extent  of  the  measurement 
volume  of  Oa  69 /mi  for  the  smaller  of  the  two  (Eq.  13),  in  reasonable  accord 
with  the  visually  estimated  value  of  ~  50  //m  using  the  cathetometer.  Tliis  value 
was  u.sed  to  calculate  the  concentration  spectrum,  at  x/d  =  100  and  He  = 

1.25  X  10"*,  plotted  as  the  solid  line  in  Fig.  3.  It  was  computed  from  the  fluorescence 
spectrum  recorded  at  the  higher  resolution,  using  Eq.  12  with  the  estimated  single¬ 
pole  transfer  function  of  Eq.  14,  at  =  r-i. 

The  effective  pole  locations  for  data  recorded  at  .r/d  =  305  were  more  difficult 
to  estimate.  At  x/d  =  305,  the  higher  relative  spatial  resolution  pushed  the  poles 
closer  to  the  noise  cross-over  point.  On  the  other  hand,  at  x/d  =  305,  the  (logarith¬ 
mic)  difference  between  the  fluorescence  and  estimated  concentration  spectra  was 
much  smaller  over  the  frequency  range  of  interest.  Figure  5  plots  the  fluorescence 
spectrum  (dotted  line)  at  x/d  =  305  and  He  =  1.2  x  10*  as  well  as  the  estimated 
concentration  spectrum  (dashed  line).  As  can  be  seen,  the  effects  of  compensation, 
in  this  case,  are  much  smaller  {cf.  difference  at,  say,  —  2.7  in  Figs.  3  and  5). 
The  estimated  concentration  spectrum  at  x/d  =  305  in  Fig.  5  is  the  one  plotted  as 
a  dashed  line  in  Fig.  3. 
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Fig.  5  Dotted  line;  Fluorescence  spectrum.  Dashed  line:  Estimated  concentration 
spectrum.  Data  recorded  at  x/d  =  305  and  J?e  =  1.2  x  10'*. 

The  jet  fluid  concentration  spectra  to  be  discussed  below  were  all  estimated 
in  this  fashion.  The  values  of  r*,  the  transfer  function  time  constant  used  in  the 
compensation  calculations,  were  fixed  for  all  the  data  measured  at  each  x/d  axial 
location  (Eq.  17,  for  measurements  at  x/d  =  100,  depending  on  which  slit  width 
was  used  to  record  the  fluorescence  data).  A  fixed  pair  of  values  was  also  used  for 
all  the  data  measured  at  x/d  =  305. 


4.  Results  and  discussion 

The  conspicuous  agreement  between  the  concentration  spectra  at  x/d  =  100 
(solid  line)  and  x/d  =  305  (dashed  line)  in  Fig.  3,  up  to  frequencies  limited  by 
signal-to-noise  ratio  considerations,  should  be  noted.  A  similar  independence  of  the 
scaled  spectra  with  downstream  location  was  also  found  to  hold  in  gas-phase  jets,^^ 
where  the  relatively  larger  diffusion  scales,  at  Sc  %  1,  and  comparable  Reynolds 
numbers  made  it  possible  to  estimate  the  concentration  spectra  with  enotigh  spatial 
resolution  directly,  obviating  the  need  for  the  compensation  scheme  employed  here. 
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In  the  case  of  the  j)resent,  liquid-phase  data,  we  should  entertain  the  possibility 
that  the  collapse  of  the  spectra  may  be  attributable  to  a  fortuitous  choice  of  the 
compensation  pole.  Some  ob.servations.  however,  are  relet-ant  here.  First,  the  close 
agreement  between  the  two  .s])ectra  also  holds  in  the  /r^  <  2  frequency  range, 
where  the  effects  of  compensation,  even  for  the  xfd  =  100  data,  are  negligible. 
As  we  will  see  below,  this  will  also  be  found  to  be  the  case  at  higher  Reynold.s 
numbers.  Second,  the  value  of  the  pole  represents  a  single  parameter.  In  contrast, 
the  collap.se  of  the  spectra  at  the  two  axial  locations,  with  substantially  different 
degrees  of  compensation  in  each  case,  is  over  5  orders  of  magnitude  in  the  power 
spectra.  Conversely,  this  collapse  supports  the  validity  of  the  single-pole  model  for 
the  leading  order  behavior  of  the  measurement  transfer  function. 


Fig.  6  Frequency-scaled  concentration  spectra.  Solid  line:  x/d  =  100  and  Re  — 
1.25  X  10"*.  Dashed  line:  x/d  =  305  and  Re  =  1.2  x  lO'*.  Dotted  line: 
reference  line  at  a  Vs  slope,  corresponding  to  a  spectrum. 


It  is  useful  to  plot  the  product  of  the  concentration  spectra  with  as  is 

commonly  done.  A  spectrum  described  by  a  -Va  power-law  yields  a  horizontal 
line  over  the  -  V.?  frequency  range  when  plotted  in  this  fashion.  The  product  of 
the  concentration  spectra  with  derived  from  the  data  at  x/d  =  100  and 
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xjd  —  305.  at  Re  =  1.25  x  10'  and  Re  =  1.2  x  10'.  arc  plotted  in  Fig.  6  a.s  the 
solid  and  dashed  line  curves,  respectively.  Also  i)lotti’d.  for  reference,  is  a  straight 
line  with  a  slope  (logarithiuic  derivative)  of  corresponding  to  the  high-Schinkit 
number  Batchelor  A-“'  theoretical  .spectrum.'' 


Fig.  7  Frequency-scaled  concentration  spectra.  Solid  lines:  xjd  =  100  data  (1.25  < 
Re  X  10“''  <  7.2).  Dashed  lines:  xjd  =  305  data  (1.2  <  Re  x  10"'  <  6.5). 
Individual  spectra  are  offset  by  — 21ogjo(f7e//7eo),  J7eo  =  1.2  x  10',  for 
clarity. 


The  solid  lines  in  Fig.  7  plot  spectra  derived  from  measurements  at  x/d  =  100. 
for  Re  X  10"'  =  1.25,  1.76,  2.55,  3.6,  5.1,  and  7.2.  The  decrease  in  the  scaled 
frequency  resolved,  in  the  x/d  =  100  spectra,  as  the  speed  of  the  flow  increased 
with  increasing  Reynolds  numlrer,  is  evident.  The  dashed  lines  plot  spectra  mea¬ 
sured  at  x/d  =  305,  for  Re  x  10"'  =  1.2.  2.4,  4.0,  and  6.5.  As  can  be  seen,  the 
highest  frequency  resolved  in  the  x/d  —  305  spectra  is  a  much  weaker  function  of 
Reynolds  number,  signal-to-noise  ratio  limitations  being  more  serious  than  spatial 
resolution  at  this  station.  Individual  specti'a  for  both  x/d  stations  are  plotted  offset 
by  — 21og,Q(Re/Reo),  with  Rcq  =  1.2  x  13',  to  aid  in  visualizing  the  evolution  of 
trends  with  Reynolds  number. 
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The  longer  characteristic  local  time  scale  at  x/d  =  305  did  not  allow  as  many 
large  scale  structures  to  be  included  in  the  record.  As  a  consc(juenc<>.  the  statistic;il 
convergence  of  the  x/d  —  305  spectra  is  not  as  good  as  for  th<'  xjd  ~  100  data. 
On  the  other  hand,  the  higher  relative  spatial  ie.solution  at  xjd  —  305  allowed  the 
spectrum  to  be  estimated  to  a  higher  (scaled)  frequeuc}-.  That  trade-off  aside,  the 
agreement  between  the  (scaled)  s])ectra  at  x/d  =  100  and  x/d  305  holds  for  all 
the  cases  for  which  data  were  recorded  at  the  same,  or  nearly  the  same.  Reynolds 
number  at  the  two  stations. 


loSloO  Tj) 

Fig.  S  Frequency-scaled  concentration  spectra  derived  from  data  at  x/d  =  100  and 
1.25  <  He  X  lO"'*  <  7.2  (no  offsets). 

Following  the  transition  out  of  the  large  scale  frequency  regime  {fTf  ~  1),  the 
spectra  appear  to  be  described  by  a  power-law  with  an  exponent  that  is  increasing 
from,  roughly,  1.2  to  1.5.  in  absolute  value,  with  increasing  Reynolds  number  (cf. 
Ref.  15,  Fig.  5.2).  This  progression  with  Reynolds  number  is  easier  to  discern  in 
Fig.  8,  which  plots  the  concentration  spectra  at  x/d  =  100,  for  He  x  10'"’  =  1.25 
(solid  line),  1.76  (dashed  line),  2.55  (dot -2-dash),  3.6  (2-dot-2-dash),  5.1  (dot 
dash),  and  7.2  (2-dot-dash),  with  no  offsets.  The  extent  of  the  power-law  regime 
can  be  seen  to  increa.se  slightly  with  increasing  Reynolds  number.  A  similarly 
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increasing  spectral  slope  (in  absolute  viilue)  with  Reynolds  number  was  also  noted 
in  measurements  of  gas-phase  spectra,'^  for  Re  x  10"'*  =  0.5,  l.C,  and  4.0. 

As  can  be  seen  by  sighting  along  the  si)ectra  in  Fig.  8,  the  power-law  region 
is  followed  by  a  different  regime  at  higher  frequencies  y<  t.  This  regime  does  not 
support  the  Batchelor  prediction"*  that  should  apply  for  over  a  decade  and  a 
half  in  frequency  in  this  case  (recall  that  Sc  ^  1.9  x  10^  here).  This  can  be  seen  in 
Fig.  6,  which  includes  a  dotted  line  with  a  reference  slope  of  -/j  .  corresponding  to  a 

spectrum.  The  spectra  at  increasing  Reynolds  number  move  even  further  fnnn 
this  slope,  as  can  be  seen  by  comparing  the  Re  =  1.2  x  10*  spectrum  with  those  at 
higher  Reynolds  numbers  in  Fig.  S.  It  should  be  noted  that  this  conclusion  extends 
to  frequencies  below  frs  ~  2.  which  are  unaffected  by  resolution  and  compensation 
considerations. 


103lO<^ 


Fig.  9  Spectrum  slope  (logarithmic  derivative).  Solid  lines:  data  at  r/d  =  100  and 
Re  =  1.25  X  10^.  Dashed  line:  x/d  =  305,  Re  =  1.2  x  10*. 


To  facilitate  the  study  of  this  higher  frequency  regime,  we  also  computed  the 
slope  of  the  spectra  (logarithmic  derivative).  A  jdot  of  spectrum  slopes,  for  the 
lower  Reynolds  number  data,  appears  in  Fig.  9.  These  were  derived  from  the  two 
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x/d  =  100.  Re  —  1.25  x  10’  spectra  iu  Fig.  3  (solid  lines)  and  the  x/d  ~  305. 
Re  =  1.2  X  10^  spectrvun  (dashed  line).  Their  coiiiparison  helps  as.sess  issues  of 
statistical  convergence  and  confidence  in  this  more  deliciite  statistic,  fus  well  a.s  the 
effects  of  spatial  resolution  and  the  applied  coinp<’n.sation. 


The  plots  in  Fig.  9  suggest  a  frcquency-dcpcndencc  of  the  sioj^e  of  the  sjtectra 
that  is  close  to  a  straight  line,  in  these  coordinates,  for  frequencies  liiglter  than 
/r^  ~  1.2.  at  this  Reynolds  number.  This  linear  behavior  extends  for  a  decade 
and  a  half  and  into  frequencies  beyond  which  the  data  tire  limited  by  resolution  and 
signal-to-noise  considerations.  The  straight  line  appears  to  be  a  good  repre.sentation 
for  frequencies  below  /r/,  2.  for  wliich  the  effects  of  comiiensation  were  negligible 

even  for  the  x/d  =  100  data  (cf.  Fig.  3).  For  frequencies  above  f  t/,  ~  2.  the  same 
straight  line  also  describes  the  behavior  for  both  tlie  x/d  =  100  and  x/d  —  305 
data,  which  were  affected  by  resolution  (and  compensation)  to  a  different  extent 
{ c/.  Figs.  3  and  5). 


Fic:.  10  Slojte  (logarithmic  derivative)  of  spectra  from  data  at  x/d  =  100  and  1.25  < 
Re  X  10“^  <  7.2.  Spectra  offset  as  in  Fig.  7.  Line  types  as  in  Fig.  8. 
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A  straight  line  for  the  logaiithmic  derivative  of  tlie  spectrum  corresponds  to  a 
spectrum  that  is  parabolic  in  log-log  coordinates,  or  log-normal  in  linear  coordinates. 
i.e., 

SAfr/^)  exp|-^  («  ln{/r;,)  +  /*]*|  .  (18) 

This  expression,  rather  than  a  power  law,  seems  to  be  the  appro[)riate  de¬ 
scription  of  our  jet  fluid  concentration  spectra  at  high  fre<picncies.  over  a  range  of 
frequencies  at  least  as  large  as  v^Sc,  i.e..  a  decade  and  a  Inilf,  in  this  case. 

Figure  10  plots  the  local  slope  (logarithmic  derivative)  of  the  spectra  in  Fig.  S. 
The  offset  scheme  employed  in  Fig.  7.  and  line-types  emploj'ed  in  Fig.  S.  were  also 
used  here.  Straight  lines  can  be  seen  to  be  a  good  representation  for  the  spectnim 
slope  at  high  frequencies  with  Reynolds  number  dependent  values  of  the  param¬ 
eters  a  and  b  in  Eq.  IS.  The  end  of  the  power-law  regime  and  the  beginning  of 
the  log-normal  range  can  be  seen  to  shift  to  higher  freciuencies  with  increasing 
Reynolds  number.  Our  data  admit  a  Kolmogorov  scaling  for  the  Reynolds  num¬ 
ber  dependeirce  of  this  transition  frequency,  i.e.,  but  corresponding  to  the 

convection  (passage  frequency)  of  a  physical  scale  roughly  SO  times  larger  than  the 
estimated  Kolmogorov  scale  (computed  using  the  estimated  mean  centerline  energy 
dissipation^^  and  the  kinematic  viscosity).  This  transition  is  not  very  well  defined, 
however,  and  other  Reynolds  number  scaling  possibilities  cannot  be  ruled  out. 


5.  Conclusions 

This  work  has  investigated  temporal  scalar  (jet  fluid  concentration)  power  spec¬ 
tra  on  the  centerline  of  high  Schmidt  number  turbulent  jets,  in  the  Reynolds  number 
range  1.2  <  Re  x  10“^*  <  7.2.  Our  .spectra  exhibit  a  power-law  regime  at  frequen¬ 
cies  above  the  local  large  scale  passage  frequency,  with  a  Reynolds  number  depen¬ 
dent  exponent  increasing  (in  absolute  value)  from,  roughly,  —1.2  to  —1.5  over  the 
Reynolds  number  range  investigated.  This  corroborates  a  similar  finding  for  gas- 
phase  jet  fluid  concentration  spectra  measured  at  comparable  Reynolds  numbers. 
At  higher  frequencies,  the  sjiectra  are  well  represented  by  a  log-normal  relation  with 
Reynolds  number  dependent  coefficients.  While  our  data  admit  a  Kolmogorov-like 
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scaling  for  the  beginning  of  the  log-nomial  region  in  the  spectrum,  i.e..  ~ 
other  possibilities  cannot  be  ruled  out. 

We  appreciate  that  our  results  are  at  odds  with  the  classical  picture  of  high 
Schmidt  number  scalar  spectra.  We  do  not  find  a  -^3  power-law  regime,  even 
though  our  measurements  were  conducted  at  Reynolds  numbers  where  such  behav¬ 
ior  has  previously  been  reported  for  high  Schmidt  number  jet  fluid  scalar  spectra.^® 
Finally,  despite  adequate  resolution  and  signal-to-noise  ratio,  our  data  do  not  sup¬ 
port  the  Batchelor  power-law  prediction.'^  Specifically,  we  found  no  constant 
A’~'  slope  at  high  frequencies  and  a  spectral  slope  that  docs  not  even  locally  attain 
a  value  of  —1. 

On  the  whole,  our  scalar  spectra  are  rather  similar  to  those  derived  by  Gargett 
from  ocean  measurements.'*  In  conjunction  with  her  data  and  analysis,  the  current 
results  raise  further  questions  about  the  universal  descriptions  of  scalar  spectra,  and 
their  applicability  to  some  of  the  canonical  flows,  such  as  turbulent  jets. 
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